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The high energy limit of QCD is controlled by the small- x part of a hadron wave- 
function. I argue that this part is universal to all hadrons and is composed of a new 
form of matter: a Colored Glass Condensate. This matter is weakly interacting at 
very small x, but is non-perturbative because of the highly occupied boson states 
which compose the condensate. Such a matter might be studied in high energy 
lepton-hadron or hadron-hadron interactions. 



1 Introduction: Parton saturation at small-x 

There has been much activity in the last few years in an attempt to understand 
the physics of nuclear and hadronic processes in the regime where Bjorken's x 
becomes very small 1-13 (and References therein). A remarkable feature about 
this regime is that the gluon density in the hadron wavefunction is so high that 
perturbation theory breaks down even for a small coupling constant, because 
of the strong non-linear effects. This is quite similar to high temperature QCD 
where the physics at the scale g 2 T is non-perijirbative because of the large 
thermal occupation numbers of the soft gluonsEJ (and Refs. therein). 

The important new phenomenon which is expected in these conditions is 
parton saturation, which is the fact that the density of partons per unit phase 
space (dN / d 2 p± dr d 2 x j_) cannot grow forever as x becomes arbitrarilly small 
(here, r = ln(l/a;)). If this growth were to occur, then the cross section for deep 
inelastic scattering at fixed Q 2 would grow unacceptably large and eventually 
violate unitarity bounds. One rather expects the cross section to approach 
some asymptotic value at high energy, that is, to saturate. It is also intuitively 
obvious that this should happen, since if the phase space density becomes too 
large, repulsive interactions are generated among the gluons, and eventually it 
will be energetically unfavorable to increase the density anyfurther. 

Saturation is expected at a phase-space density of order l/a s , a value 
typical of condensates trll This, together with the fact that gluons are massless 
bosons, leads naturally to the expectation that the saturated gluons form a new 
form of matter, which is a Bose condensate. Since gluons carry color which is 
a gauge-dependent quantity, any gauge-invariant formulation will neccessarily 
involve an average over all colors, to restore the invariance. This averaging 
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procedure bears a focmal resemblance to the averaging over background fields 
done for spip_glassesE3. The new matter is therefore called the Colored Glass 
CondensateEl(CGC). 

In what follows, I will briefly explain the assumptions which justify this 
physical picture, and the mathematical formalism which is used to describe it. 

2 The Colored Glass Condensate 

The CGC picture holds, strictly speaking, only in the infinite-momentum 
frame, where the hadron propagates almost at the speed of light, and, by 
Lorentz contraction, appears to the external probe as an infinitesimally thin 
two-dimensional sheet located at z = t, or x~ = 0. In what follows, I shall 
use light-cone (LC) vector notations: for an arbitrary 4- vector v^, I write 
v» = (v+,v~,v±), with v+ = (l/\/2)(t; +v 3 ), v- = (1/V2)0° - v 3 ), and 
vj_ = (v ,v ). The dot product reads then: p ■ x = p + x~ + p~x + — p± ■ x±; 
p~ and p + are, respectively, the LC energy and longitudinal momentum; cor- 
respondingly, x + and x~ are the LC time and longitudinal coordinate. 

In the infinite-momentum frame, the parton interpretation makes sense 
and deep inelastic scattering (DIS) proceeds via the instantaneous absorbtion 
of the virtual photon 7* (with 4-momentum q M ) by some parton in the hadron. 
The Bjorken xb variable is defined as xb = Q 2 /2P-q, where Q 2 = —q^q^, and 
pp- — S^ + P + with large P + is the hadron 4-momentum. By kinematics, xb 
coincides with the longitudinal momentum fraction x = p + / P + of the struck 
parton: ig = x. 

At x <C 1, the gluon density increases faster, and is the driving force 
towards saturation. The dynamics responsible for this increase is the BFKL 
evolution til, as I briefly recall now (see also Ref. tl and Rcfs. therein): 

Let me call soft, respectively fast, a parton having small, respectively, 
large longitudinal momentum (this separation will be made more precise in 
a moment). A soft gluon, with p + = xP + <C P + , is a shortlived excitation 
which is typically radiated by a fast parton (e.g., a valence quark) with a larger 
longitudinal momentum k + p + , and thus a longer lifetime. Indeed, by the 
uncertainty principle, the lifetime of the parton system in Fig. l.a is: 

. 1 I 1 2p+ 

Ax = ~ — = — — oc x, (1) 

Efe-p + £p - Sk S p pj_ 

where = k 2 _/2k + is the LC energy of the on-shell excitation with momen- 
tum k = (fc + ,kjj, and I have used the fact that, for comparable transverse 
momenta k± and p±, e p S> Ek,£k-p- To conclude, softer partons have larger 
energies, and therefore shorter lifetimes. 
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Figure 1: a) Soft gluon emission by a fast parton; b) a second emission; c) a gluon cascade. 

The lowest-order process in Fig. l.a is amended by radiative corrections 
enhanced by the large rapidity gap At = ln(fc + /p + ) ~ ln(l/x). (I assume here 
that ln(l/x) ^> 1.) For instance, the probability for the emission of a second 
gluon with momentum in the range k + > > p + is (cf. Fig. l.b) : 



and becomes exceedingly large as x — > 0. It is then highly probable that more 
gluons will be emitted along the way, thus giving birth to the gluon cascade 
depicted in Fig. I.e. For a fixed number N of gluons in this cascade, the largest 
contribution, of order (a s ln(l/a;)) , comes from the kinematical domain where 



Other momentum ordcrings give contributions which are suppressed by, at 
least, one factor of l/ln(l/a;), and thus can be neglected to leading logarith- 
mic accuracy (LLA). For the dominant contribution in eq. (|J), the number of 
radiated gluons increases exponentially with At: N(x) ~ exp{Aa s ln(l/a;)}, 
with constant A. This is a coherency effect, consequence of the separation of 
scales in eq. (|J), and of the bosonic nature of the gluons: 

Indeed, because of its short lifetime, the soft gluon at the lower end of 
the cascade "sees" the N previous gluons as a frozen color charge distribution, 




k + = fc+ > fc+ > fc+ > ■ • ■ > fc+ = p + . 



(3) 
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with an average color charge Q = \J (Q a Qa) oc N. Thus, the iVth gluon is 
emitted coherently off the color charge fluctuations of the previously emitted 
gluons, with a differential probability (compare to eq. (j^)) : 

dP N oc a s N dr N , (4) 

which implies that N(t) ~ e AasT , as anticipated. Then, also the gluon density 
xG(x, Q 2 ) = idN/dr) oc e AotBT grows exponentially^ with r = \n(l/x). 

Thus, the BFKL picture is that of an unstable growth of the color charge 
fluctuations as x becomes smaller and smaller. However, this evolution assumes 
the radiated gluons to behave as free particles, so it ceases to be valid at very 
low x, where the gluon density becomes so large that the radiated gluons 
overlap each other in the transverse plane and start interacting. This is the 
onset of saturation. 

This is also the regime where the description in terms of a Colored Glass 
Condensate becomes appropriate: Because of the hierarchy of scales in eq. (||), 
the soft gluons "see" the fast partons as an effective color charge which is static 
(i.e., independent of x + ), and localized near the LC (i.e., at x~ =0), with a 
random density p a (x~ , x±). (This is random since the soft gluons can belong to 
different cascades, and the instantaneous configuration of the cascades inside 
the hadron is random. Of course, the integrated color charge over the whole 
hadron must vanish, because of confinement.) The spatial correlations of the 
effective charge p a {x~~ , x±) reflect the quantum dynamics of the fast partons, 
and are encoded in a statistical weight function W[/o]. 

This color charge acts as a source for the soft gluons which, because of their 
large occupation numbers, can be treated in the classical approximation. This 
suggests the following effective theory for_the soft gluon correlators, originally 
proposed by McLerran and Venugopalancl: First, one solves the classical Yang- 
Mills equations with the source p a , and in the LC-gauge^ A + — : 

[D V ,F"»] = 6»+p a (x-,x ± ). (5) 

In the saturation regime, A ~ p ~ \Jg (see Sect. 3.2 below), and eq. (||) must 
be solved exactly: the classical problem is fully non-perturbative. The non- 
linear effects correspond to interactions among the soft gluons, to all orders. 
The corresponding solution A 1 - 1 = A^[p] will be constructed in Sect. 3.1. Then, 
the soft correlation functions of interest are obtained as: 

(AZ(x)A v b (y) ■ ■ -)a = J Vp W A [p] A%(x)Al(y) (6) 

C A more refined treatment, using the BFKL equation Q, gives A = 4N C \n2/ir. 
d This gauge choice will be motivated in Sect. 3.2. 
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Note the scale A = A + in eq. (^): this is the separation scale between fast 
(p + > A+) and soft (p + < A+) degrees of freedom. Clearly, the structure and 
the correlations of p will depend upon A (e.g., p has support at x~ < 1/A+, 
and W[p] = Wa[p])- Moreover, the soft correlations must be evaluated at lon- 
gitudinal momenta p + which are not too small as compared to A + : otherwise, 
the quantum corrections due to the "semi-fast" gluons with k + momenta in 
the range p + < k + < A + would be relatively large, of 0(a s ln(A + /p + )), since 
enhanced by the large rapidity interval At = ln(A + /p + ) ^> 1. 

Remarkably, the equations (§|)-(§|) are those for a glass (here, a colored 
glass): There is a source, and the source is averaged over. This is entirely 
analogous to what is done for spin glasses when averaging over stochastic 
magnetic fields E£l. By computing the two-point function as above, one has 
foundO'EI a saturation regime where {A l a A l a )A oc l/a s , a value typical for a 
condensate (see also Sect. 3.2 below). We thus conclude that the matter 
which describes the small x part of a hadron wavefunction is a Colored Glass 
Condensate. 

3 Saturation in the McLerran-Venugopalan model 

In this section, I present the solution to eq. (||) and explore its consequences 
for the saturation of gluons. 

3.1 The non-Abelian Weizsacker-Williams field 

Consider first the Abelian version of eq. (|^), as a warming up: 

d v F v » = 5» + p(x), (7) 

where A+ = 0, = d^A" - d u 'A" , and x= (x _ ,xx). Eq. (0) can be easily 
solved in momentum space, with the result that A~ = and 

A*{p) = -^e^pA. (8) 
p+ P z ± 

Here, one needs a prescription to handle the pole at p + = 0: different prescrip- 
tions correspond to different boundary conditions in x~ . For instance, with 
the "retarded" prescription^] 1/ (p + +ie), one obtains 

A i (x~,x±)= [ dy- d l a{y-,x ± ), (9) 



e Note that the "retardation" refers here to x , and not to the LC time av 
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with a satisfying — V^_a;(af) = p(x). 

The vector potential (||) is static, d~A l = (dA l /dx + ) = 0, and defines a 
two-dimensional pure gauge: T 1 ^ — 0. But, of course, this is not a pure gauge 
in four dimensions, since the electric field T l+ = — d + A l = d l a is non-zero. 
The above solution, also known as the Weizsacker- Williams field, is the analog 
of the Coulomb field in the infinite momentum frame. 

It turns out that the corresponding non-Abelian solution has a similar 
structure: Specifically, the solution A 1 to eq. (||) in the LC gauge (A + = 0) can 
be chosen such as A~ — 0, and A is static and a two-dimensional pure gauge: 
jnj = o. This can be written as follows (for retarded boundary conditions): 

A^x-.xx) = (i/g)U{x-,x ± )d l U' l {x-,x ± ) (10) 

U^(x~, xj_) = P cxp <|zg J dz~ a(z~,x±)^j , (11) 

- V\a{x) = p{x) = U\x)p{x) U(x), (12) 

where P is the usual path-ordering operator, and p(x) is the value of the 
classical source in the covariant gauge d^A^ = 0. 

Note that, while a is linearly related to p, its relation to the LC-gauge 
source p is more complicated, since the gauge rotations U and W are them- 
selves functionals of a, cf. eq. (^). Thus, in order to perform the average in 
eq. (|^) , it is more convenient to use p (rather than p) as the independent vari- 
able. This is possible because the measure and the weight function in eq. (^|) 
are gauge invariant; e.g., Wa[p] = \p) . Then, the non-Abelian solution 
A l [p\ is known explicitly, via eqs. and eq. (|) is replaced by 

(A\x + ,x)Ai{x+,y) ■ ■ -)a = J Vp W A [p] A s [p] Al[p] (13) 

In particular, this equation shows that only the equal-time (x + — y + ) corre- 
lators of the soft transverse (jjl = i, i = 1,2) fields can be computed in this 
model; but these are precisely the gluon correlators which are probed in DIS. 

Recall finally that the source p is created by the fast partons with \p + \ > 
A + , so its support is restricted to ^ 1/A + . The small- a; external probe, 
on the other hand, is sensitive only to the gross features of the color fields A 1 
over large distances \x~\ ^> 1/A+, where one can replace eq. ( |TT| ) with 

U\x~ , x±_) » #(x~)Pexp < ig / dz~ a(z~ , x±) > = 0(x~) v\xx)^ (14) 
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and therefore (cf. eq. (pX|)): 

A 1 ^6{ x -)-v{d l v ] ) ee 9(x-)Ai (x 1 _), F+ ^ -5(x-)A i x (x±_), (15) 
9 

where, strictly speaking, the 5-function is localized at x~ < 1/A+. 
3.2 Gluon distribution function and saturation 

Consider the simple approximation where, the weight function for p, which is 
not yet known, is taken to be a GaussianEl : 

^A[/5]^exp|-i y^Hg}, ( 16 ) 

where £^ is the total color charge density squared of the partons with p + > 
A + . By using this approximation, it has been possible to compute the gluon 
distribution function xG(x, Q 2 ) in the MV model EH3, as I recall now: 

By definition, G{x, Q 2 ) is the number of gluons in the hadron wavefunction 
having longitudinal momentum k + = xP + , and transverse momentum less 
than Q. In the LC gauge, this is simply related to the Fock-spafie gluon 
distribution function^], and therefore to the gluon two-point functionlil: 

x (^A a {x + ,k+^ A _)A a {x + ,-k + ,-^ A _)), (17) 

where the brackets denote the average over the hadron wavefunction. In the 
LC-gauge gauge, F l+ (k) = ik + A l (k), and therefore (with k + = xP + ) : 

xG(x,Q 2 ) = i| |^e(Q 2 -kl)(F?(x + ,k)Fi + (x + ,-k)). (18) 

If the integration over k± were unrestricted, this quantity would be manifestly 
gauge- invariant. But even for a finite Q 2 , this has a gauge- invariant meaning 
when evaluated on the non-Abelian Weizsacker- Williams field in Sect. 3.1LL3. 
In this approximation, F l+ (x + ,k) w T 1+ {k±) — — A l oa (k±) (cf. eq. (jlq)), and 

xG(x,Q 2 )*R 2 J Q ^ J d 2 x x e- ik ^(A£(0)A™( Xl _)) A , (19) 

f Of course, the Fock-space gluon distribution can be defined in any gauge; biitipfeiis only in 
the LC-gauge that the definition (|17j) can be given a gauge invariant meaningliltj. 
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where R is the hadron radius (I have assumed homogeneity in the transverse 
plane, for simplicity), and the average is to be understood in the sense of 
eq. ( |l3| ) with A + = xP + . Thus, the r.h.s. of eq. ( |l9| ) is still dependent on x, 
but only via the respective dependence of the weight function for p. 

With the Gaussian weight function ( |l6|) , and the non- linear classicaLsolu- 
tion in Sect. 3.1, the gluon distribution (|19|) can be computed exactly □'□: 



N 2 



1 1 



1TOi s N c 



-x 2 ± ln(x 2 ± A% GD )Q 2 s /4 



(20) 



where N c is the number of colors, and Q s oc a s ^\ is the saturation momen- 
tum and is a function of A + , that is, of Bjorken's x. This equation displays 
saturation: the vector potential never becomes larger than A 1 ~ 1/g. This 
is the maximal occupation number permitted for a classical field, since larger 
occupation numbers are blocked by repulsive interactions of the gluon field. 

This interpretation can be made sharper by going to momentum space: If 
N(k±) is the Fourier transform of (20) [this is the same as (dN/d 2 k± dr d 2 x±), 
the gluon density per unit rapidity and unit transverse phase-space], then: 



N(k ± ) oc a s (Q 2 Jkl) for k\ > Q 2 S 



(21) 



which is the normal perturbative behavior, but 



N(k A 



1 



In 



for 



A 2 

1V QCD 



< fci < Ql, 



(22) 



which shows a much slower increase, i.e., saturation, at low momenta. 

Note, however, that the above argument is not rigurous, since the local 
Gaussian form for W[p] in eq. (16) is valid only at sufficiently large transverse 
momentum scales so that the effects of high gluon density are small. It is 
therefore important to verify if saturation comes up similarly with a more 
realistic form for the weight function, as obtained after including the quantum 
evolution in x. This will be discussed now. 



4 The non-linear evolution equation 

Because the separation of scales is only logarithmic, the effective theory (||)-(||) 
with scale A + applies only to gluon correlations at a scale p + slightly below 
A + . If one is interested in correlations at the softer scale bA + with !) « 1, 

9 In relativistic heavy ion collisions, one expectsE) Q a ~ 1 GeV atRHIC, and Q s ~ 2 — 3 GeV 
at LHC. Thus, at least at LHC, the kinematical window in eq. (E3) should be non-negligible. 
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then, to LLA, one has to include also the corrections of order a s ln(l/6) due 
to the semi-fast quantum fluctuations with longitudinal momenta in the strip 



bA+ < \k+\ < A+. (23) 

Together with the dynamical information already contained in the effective 
theory at scale A + ("Theory 1"), these additional corrections will determine 
the effective theory at the softer scale 6A + ("Theory 11"). This suggests an 
iterative construction of the effective theory where the quantum fluctuations 
are integrated out in layers of k + , down to the physical scale of interest. At 
each step in this procedure, the quantum corrections must be computed to 
leading order in a s ln(I/6) (LLA), but to all orders in the strong background 
field A 1 ~ 1/g produced by the color source p at the previous step. 

To compute quantum correctioaSj-oae needs the quantum generalization of 
the McLerran-Venugopalan modelBBli3. This is obtained by replacing eq. (||) 
with (below, the gauge condition A + = is implicit): 

J { j A T>A e %s l A > A J 

which involves two functional integrals: a quantum path integral over the 
soft (k + < A + ) gluon fields A 11 , which defines quantum expectation values 
at fixed p, and a classical average over p, with weight function Wa[p]- (This 
double averaging is similar to the one performed for.-spin systems at finite 
temperature and in a random external magnetic field ta.) Unlike eq. @, the 
2-point function given by eq. (^4|) is independent of the arbitrary separation 
scale A + : the cutoff dependence of the quantum loops cancels against the 
corresponding dependence of the classical weight function Wa [p] ■ 

The action S[A, p] is chosen such as to be gauge-invariant and repro- 
duce the classical equations of motion (|^) in the saddle point approximation 
SS/SA^ = 0. ft reads!: S = S YM + S w , where S YM = J ^(-F^JA) is the 
Yang-Mills action, and Sw is a gauge-invariant generalization of the cikonal 
vertex J d xpA~ (with T denoting time-ordering of color matrices) : 

S w = J d 3 xTr {p(x)Texp(ig J dx + A-(x + ,x^j | . (25) 

Theories I and 11 are defined as in eq. (|24|), but with separation scales 
equal to A + , and 6A + , respectively. The difference AW = WbA — Wa can be 
obtained by matching calculations of gluon correlations at the scale p + bA + 
in the two theoricsta. In Theory II, and to lowest order in a s , these correlations 
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are found at tree level, i.e., in the classical approximation. In Theory I, and 
to the same accuracy, they involve also the logarithmically enhanced quantum 
corrections due to the "semi- fast" fluctuations defined by eq. ((2^). The result 
can be expressed as an evolution equation)^] for W T [p] (with r = ln(P + /A + )) 
with respect to variations in r, originally derived in Ref. □. This reads: 

dWM fl ^ [W TXxy ]-^[W T a x] \. (26) 



<9t [2 Sp x Sp y S p. 

in condensed notations where (with x ~ = 1/A + = c t 2;q~), e.g., 

t— [W T o x ] = / d 2 x ± —— W T a a {x ± ) , (27) 

op x J op a (x T , x± ) L J 

and the functions a a (x±) and Xab(^±,y±) are the one- and two-point correla- 
tors of the color charge dp a (x) of the semi-fast quantum fluctuations: 

a s ln~<7 a (xj_) = J dx~ {5p a (x + , x~ , x±_)) p , (28) 

a s \n^-Xab(xA_,y±_) = J dx~ J dy~ {6p a (x + , x~ , x_i_)5p b (x + , y~ , y±)) p . 

Thus, a and \ are (generally non-linear) functionals of p obtained by averag- 
ing over the semi-fast quantum fluctuations in the presence of the tree-level 
fields and sources, A 1 and p. To the order of interest, the semi- fast fields can 
be integrated out in the Gaussian approximation; thus, the Feynman graphs 
contributing to a and x involve, at most, one loop. Some typical contributions 
to lowest order in p are shown in Fig. |[ As one can see on these figures, a 
is a one-loop virtual correction to the tree-level source and describes vacuum 
polarisation, while x is a tree diagram which describes the real emission of a 
semi-fast gluon (cf. Fig l.b). 

The functional equation ( p6[ ) is equivalent to an infinite hierarchy of or- 
dinary equations for the correlators of the charge density. For instance, by 
multiplying eq. ( p(| ) with p x p y and functionally integrating over p, one obtains 
an evolution equation for the two-point function: 

T (PxPv)t = Ol s (x X y + pxCTy + (T x py)r , (29) 

ctr 

which in general involves also the higher n-point functions, via a and x- But 
in the limit of a weak source, i.e., with a and x computed to lowest order 

h This equation insures that correlation functions like ( p^j ) are independent of A+. 
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Figure 2: Contributions to a (a) and \ ( D ) to lowest order in p. Wavy lines with a blob 
denote insertions of p\ the continuous lines are propagators of the semi-fast gluons; the 
external, dotted, lines carry soft momenta. 



in p. this, becomes a closed equation for (pp) which coincides with the BFKL 
equations, as necessary on physical grounds. This is a very non-trivial check of 
the effective theory in eqs. (|2^)-(p5|). Specifically, the linear (in p) contribution 
to a depicted in Fig. 2. a is evaluated as: 



(30) 



with the "form factor" 



p —ibA + x p —iA Jr x 

F(x~) = 6{x-) - , (31) 

x 

which has support at 1/A+ ^ x~ ^ 1/6A+, and yields the expected logarithmic 
enhancement (cf. eq. (|2^)) after the integration over x~ : 

dx-F(x-) = ln T . (32) 
b 



The kernel of the two-dimensional integral in eq. (30) can be recognized as the 
virtual part of the BFKL kernel; the corresponding^ real part is generated by x 
when evaluated to lowest order in p (cf. Fig. 2.b)tl 

The general, non-linear, expressions_for a and \ valkl in the saturation 
regime will be presented somewhere elseEj. (See also Ref. Ej for an alternative 
calculation, with different results though.) Here, I would like only to empha- 
size the longitudinal structure of the quantum corrections, which is already 
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manifest on the linear approximation (|3(])-(|3l]) : Since generated by modes 
with k + momenta in the strip (p3[), the induced source (Sp) is localized at 
1/A+ < x~ < 1/6A+, that is, on top of the tree- level source p (which has 
support at ^ x~ 1/A+). Because of that, the functional derivatives in 
eqs. (p6|)-(p7|) are to be evaluated at x~ = x~ = 1/A+ : the evolution from 
W T [p] to Wr+drip] is due to changes in p within the rapidity interval (r, r+dr). 

In fact, the precise longitudinal picture depends upon the ie prescription 
for the "axial" pole 1 /p + in the LC-gauge_propagator of the semi- fast gluons^]. 
The result in eq. @ has been obtained E3 by using the retarded prescription 
+ is), for consistency with the classical solution (jl0|)-(|l2|). If the ad- 
vanced prescription l/(p + — is) was used instead, the corresponding (5p a (x)) 
would be located at negative x~ , but its integral over x~ — which defines 
a(x±), cf. eq. (|28| ) — would nevertheless be the same. That is, the BFKL 
equation is obtained independently of the axial prescription. 

This is, however, specific to the linearized (or BFKL) limit of the evolution 
equation. In general, the non-linear effects depend upon the gauge condition- as 
shown by the explicit-calculations of a and x using different prescriptionsEJilj. 
The simplest results E3 are obtained by using the retarded, or the advanced, 
prescriptions alluded to before. With a retarded prescription, both the classical 
mean field A a , eq. (|l5|), and the induced source (Sp a ), eq. (|30|), sit at x~ > 0, 
or z < t ; that is, the soft gluons fields are behind their source, the hadron 
(which is located at z = t). It is intuitively plausible that within this picture 
there are no initial state gluon interactions. This is the counterpart of the 
conclusion of Mueller and KovchegovQfl who used an advanced prescription 
l/(p + — ie), and showed that all the final state gluon interactions disappeared. 

These two prescriptions are simple in that one can either put gluon interac- 
tions in the final or initial state. Other gauge conditions such as Leibbrandt- 
Mandelstam or principle value do not have this simple feature, and lead to 
other complications for computations as well Ej. Of course, gauge invariant 
quantities must come out the same whatever gauge is used for their compu- 
tation. But the classical color charge p a (x) is not a physical observable by 
itself, but just a convenient tool to summarize the (generally, gauge- variant) 
correlations inherited by the soft gluons from the fast-oaiytons. The discrep- 
ancies between the results for a and x found in Refs. EJ'EJfij and Ref. 113 may 
be ultimately attributed to the different methods used to fix the gauge. In 
fact, our results £3 are consistent with some previous results by Balitskiitj and 
KovchegovQ, to which they reduce in the large N c limit. 

To conclude, the non-linear evolution eouation ( p6| ) together with the ex- 
pressions for a and x to be presented in Ref. Il3 provide a theoretical framework 

! This prescription is necessary for fixing the residual gauge freedom in the LC-gauge. 
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to perform calculations in the saturation regime. Further studies within this 
framework should lead to a quantitative understanding of unitarity in high 
energjiL scattering, and of the structure of final states and multiparticle produc- 
tionEm. This would have important applications to the theory of deep inelastic 
scattering, hadronic interactions and ultrarelativistic nuclear collisions. 
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